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Abstract 

We discuss the entropy of a circular polymer under a topological constraint. We call it the topological 
entropy of the polymer, in short. A ring polymer does not change its topology (knot type) under any 
thermal fluctuations. Through numerical simulations using some knot invariants, we show that the 
topological entropy of a stiff ring polymer with a fixed knot is described by a scaling formula as a 
function of the thickness and length of the circular chain. The result is consistent with the viewpoint 
that for stiff polymers such as DNAs, the length and diameter of the chains should play a central role in 
their statistical and dynamical properties. Furthermore, we show that the new formula extends a known 
theoretical formula for DNA knots. 

keyword topological entropy, circular polymers, random knotting, self-avoiding polygons, knot invari- 
ants, DNA knots 



1 Introduction 

In the last fifteen years, various knotted ring- 
polymers are synthesized and observed in experi- 
ments of chemistry and biology. |l| |[ ^ Once a 
ring polymer is formed, its topological state, which 
is given by a knot, is unique and invariant. One of 
unsolved problems in statistical physics of macro- 
molecules is how to formulate the topological con- 
straint on ring polymers theoretically so that we 
can investigate topological effects on their statis- 
tical and dynamical properties. This problem was 
addressed by Delbriick in the 60's @, |5). Since then, 
several numerical simulations have been performed 
|,g,|,g|ol0,0||0[l||l6l[l|. The topo- 
logical constraint on a ring-polymer is nontrivial. It 
may severely restrict the available degrees of free- 
dom in the configuration space, and then lead to a 
large reduction on its entropy. 

Let us formulate the topological problem, explic- 
itly. We consider a ring polymer in good solution, 
whose degree of polymerization corresponds to N 
units of the Kuhn statistical length [Q. Then, a 
spatial configuration of the ring polymer can be 
approximated by a self-avoiding polygon with N 



polygonal nodes (or N bonds, N vertices). Here, 
the length of polygonal segments is given by the 
Kuhn length. 

Let us assume the number Wk{N) of all possi- 
ble A^-noded self-avoiding polygons with a knot K. 
Then, the topological entropy Sk{N) for the knot 
K is given by Sk{N) = fee ^ogWxiN). Here ks is 
the Boltzmann constant. We remark that the en- 
tropy 5'(A^) without any topological constraint is 
given by S{N) = \ogW{N), where W{N) is the 
number of all self-avoiding polygons with A'' nodes: 
W{N) = J2k^k{N). Let us introduce the knot- 
ting probability Pk{N) of a knot K. We define it 
by the probability that a given configuration of a 
self-avoiding polygon with A^ nodes is equivalent to 
the knot K. It is clear that the knotting proba- 
bility is given by Pk[N) = Wk{N)/W{N). Thus, 
the decrease of the entropy of the ring polymer due 
to the topological constraint that it should keep its 
knot type K is expressed in terms of the knotting 
probability for the knot K: 

S{N) ~ Sk{N) = kB \og W{N) - fee log Wk{N) 
= fcBlog(l/FA'(A^)). (1.1) 
Recently, randomly knotted DNA rings are syn- 
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thesized in experiments, and the fractions of some 
knotted species are measured through the elec- 
trophoretic separation. Here we remark 

that the fraction of DNAs with a knot K corre- 
sponds to the knotting probabihty of the knot K. 
The measurement of the fractions of knotted DNAs 
is particularly interesting because it gives an ex- 
cellent method for determining the effective diame- 
ter of DNAs, which are negatively charged and sur- 
rounded by the clouds of counter ions. The effec- 
tive diameter of DNAs in an electrolyte solution de- 
pends strongly on the concentration of the counter 
ions. However, it is nontrivial to calculate the ef- 
fective diameter of DNAs by taking into account 
the electric double layers around the cylindrical seg- 
ments. |2^, Furthermore, DNA chains are rather 
stiff; the persistent length (i.e., the Kuhn statistical 
length) of DNA chains is given by 50 ~ 100 nm, 
depending on the ionic concentrations [2^, which 
is much larger than that of polystylenes. [|l8[ 

Thermodynamical properties of DNA chains in 
electrolyte solutions are studied by introducing 
some models of cylindrical self-avoiding walks. For 
DNA rings, we introduce a model of self-avoiding 
polygons consisting of iV cylindrical segments with 
cylindrical radius r where their length is given by 
the Kuhn length h. Then, let us consider the prob- 
ability Pknotted{N,r) of a cylindrical polygon be- 
ing knotted (i.e., being a non-trivial knot). For 
the algorithms of cylindrical self-avoiding polygons 
it was studied through simulations how the 
probability PknottediN , r) should depend on the ra- 
dius r. The numerical results were analyzed by as- 
suming an exponential dependence in the following: 

PknottediN, r) = PknottediN, 0)exp(^-~-f^^ . (1.2) 

Here 7 is a constant to be determined by the simu- 
lations. Thus, by measuring the fraction of knotted 



DNAs and applying the theoretical formula (1.2) to 
them, the effective diameters of DNA rings in some 
electrolyte solutions are determined. |l^, |2^ 

In this paper, we introduce a theoretical formula 
describing the knotting probability of a given knot 
for stiff ring polymers such as circular DNAs. It de- 
scribes how the knotting probability depends on the 
length and radius of the ring polymers. We formu- 
late an algorithm of self-avoiding polygons consist- 
ing of cylindrical segments with radius r and unit 
length. The algorithm is based on the dimer- 
ization algorithm, and we call the algorithm the 
cylindrical ring-dimerization method, or the dimer- 
ization method, for short. The method should be 
useful for studying stiff ring polymers. In fact, it is 
closely related to the wormlike chain model for poly- 
electrolytes in electrolyte solutions. Through 
numerical simulations with the method, we eval- 
uate the knotting probabilities for the cylindrical 
self-avoiding polygons with several different num- 
bers N of polygonal nodes and different values of 



the cylindrical radius r. Then, we show that the 
new formula gives good fitting curves to the esti- 
mates of the knotting probabilities. In fact, the 
fitting curves fit to the data for the cylindrical poly- 
gons of very small as well as very large (even asymp- 
totically large) numbers N of nodes, and of all the 
different values of the radius r. Furthermore, we 
show that it extends the known formula ( |1.2| ) for 
the DNA knots. Finally, we discuss some conse- 
quences of the theoretical formula of the knotting 
probability, which may be useful for future study of 
knotted circular DNAs. 

The model of cylindrical self-avoiding polygons 
should be valid for any real ring polymer in good 
solution. According to the standard two-parameter 
theory of polymers, it is expected that any statis- 
tical property of polymers can be described by the 
length and thickness of the chains. |l^ The radius 
of cylindrical segments corresponds to the thick- 
ness of ring polymers, which expresses the excluded- 
volume effect. In appendix A, some algorithms of 
the cylindrical self-avoiding polygons are shown. 

The paper is organized as follows: In §2 we ex- 
plain the cylindrical ring-dimerization method of 
our numerical experiment of random knotting. In 
§3, employing the algorithm discussed in §2, we 
produce a large number of self-avoiding polygons 
consisting of cylinders, and present numerical esti- 
mates of the knotting probabilities of several dif- 
ferent knots such as the trivial knot, non-trivial 
prime knots (3i, 4i, 5i, 82), and nontrivial compos- 
ite knots (3ipi, 3ipitt3i). Then, we find that they 
are well described by the new theoretical formula 
of the knotting probability of the cylindrical self- 
avoiding polygons constructed by the cylindrical 
ring-dimerization method. In §4, we show that the 



new formula generalizes the known formula (1.2) 
Finally, we give some concluding remarks in §5. 



2 Methods of Numerical Sim- 
ulations 

2.1 Overview of the numerical ex- 
periment of random knotting 

We first discuss the outline of our numerical experi- 
ment. The knotting probability is evaluated by the 
following processes: (1) We construct a large num- 
ber, say M, of self-avoiding polygons of N cylin- 
drical segments with radius r and unit length; (2) 
Projecting the three-dimensional configurations of 
the self-avoiding polygons onto a plane, we make 
their knot diagrams; (3) Calculating some knot in- 
variants for the knot diagrams, we enumerate the 
number Mk of such polygons that have the same 
set of values of the knot invariants with a knot K. 

For a given number N of polygonal nodes, we pro- 
duce M polygons with the algorithm of the dimer- 
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ization method of self-avoiding polygons. In this 
paper, we set M — 10'*. We recall that Mk denotes 
the number of such polygons that have the same 
set of values of the knot invariants with the knot 
K. Then, it is clear that the knotting probability 
Pk{N) is given by the following 

Pk{N) = ^ (2.1) 

We note that the algorithm of the dimerization 
method will be explained in §2.3 and §2.4 (See also 
§A.l). 

2.2 Knot invariants for detecting 
knot types 

In the numerical simulations of the paper, the tool 
for detecting the knot type of a given polygon is a 
set of two knot invariants in the following: the de- 
termininant A/<(— 1) of knot and the Vassiliev-type 
invariant V2{K) of the second degree. The values 
of the invariants for some typical knots are given in 
Table 0. It is remarked that the two knot invariants 
have definite advantages for practical purposes. In 
fact, there exist some algorithms by which we can 
calculate the two knot invariants in polynomial time 
with respect to the number of crossing points of a 
given knot diagram. ||2^, |2^ 

2.3 Cylindrical ring-dimerization 
method 

Let us discuss the method |2^ for constructing self- 
avoiding polygons consisting of cylinders with ra- 
dius r of unit length. It is based on the dimerization 
algorithm. We recall that it is called the cylin- 
drical ring-dimerization method, or the dimeriza- 
tion method, for short. The method consists of the 
following processes: (1) We generate a set of chains 
with cylindrical segments by the dimerization al- 
gorithm; (2) We construct polygons by connecting 
two cylindrical self-avoiding chains with the method 
of Ref. ||l^ , where we also calculate the statistical 
weight related to the probability of successful con- 
catenation. 

Let us explain the dimerization algorithms for 
self-avoiding walks, and then for self-avoiding poly- 
gons. Under the dimerization algorithm P7| , we 
make a chain by connecting two given subchains 
if they have no "overlap" ; if they have an "over- 
lap" , then we throw away both of the two subchains 
and we take a new pair of subchains to try. In the 
method [ p^ for making self-avoiding polygons, first 
we pick up a compatible pair of self-avoiding walks, 
and then we check whether they have any "overlap" 
or not. If not, we make a new ring by connecting 
the ends of the two chains; if there is an "overlap" , 
then we throw them away and take a new pair of 



chains. In the concatenating method, we also calcu- 
late the statistical weight related to the probability 
of successful concatenation. 

One of the key points of our algorithm is the con- 
dition when given two cylindrical segments have an 
"overlap". We define the condition of " overlap" 
between a given pair of segments as follows. First, 
there is no "overlap" between any pair of adjacent 
segments, i.e., we neglect the thickness of segments 
for any pair of adjacent cylindrical segments. Thus, 
there is no constraint between any next-neighboring 
segments in our model of cylindrical self-avoiding 
polygons. Second, two segments have no "overlap" 
if the minimum distance between the central line 
segments of cylinders for any given pair of unad- 
jacent cylinders is larger than the diameter 2r of 
cylinders. Here, we have defined the central line 
segment of a cylinder by the line segment between 
the centers of the upper and lower disks of the cylin- 
der. For an illustration, we give in Appendix B the 
algorithm for checking whether any given two cylin- 
drical segments have an "overlap" or not. 

The cylinder condition is important in practical 
applications. We recall that DNA chain as a poly- 
electrolyte in an electrolyte solution can be modeled 
by the model of cylindrical self-avoiding polygons 
since the negatively charged DNA segments with 
counter ions shielded around them can be approx- 
imated effectively as impermeable cylinders of the 
Kuhn statistical length with the radius given by the 
Debye screening length. 

2.4 Some details of the dimerization 
algorithm 

Let us describe the algorithm of dimerization, more 
precisely, p?! First, we construct a very large num- 
ber, say Mq, of short self-avoiding chains with 
cylindrical segments, randomly by a direct method. 
Then, we pick up a pair of short chains out of the 
Mq generated chains randomly, and see whether 
there is any overlap or not between the chains; we 
investigate all unadjacent pairs of cylindrical seg- 
ments of the chains. If there is an overlap, then 
we give up the chains and consider a new pair of 
short chains from the beginning. If there is no over- 
lap, then we make a longer chain with length 2A'o 
by connecting the end of one chain to the head of 
another one. After constructing Mo/2 chains with 
length 2Nq by the method, we repeat this construc- 
tion for the Mo/2 longer chains. The scheme of the 
dimerization algorithm is shown in Fig. 3, which is 
given by a binary tree. Thus, we can construct very 
long self-avoiding polygons with the cylinder radius 
r, quite efficiently. In appendix A, we explain all 
the detail of the algorithm of the cylindrical ring- 
dimerization, explicitly. 
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Knotting Probability 
Its Fitting Formula 



and 



3.1 A brief review on knotting prob- 
abilities of ring polymers 

We discuss some important known results on the 
statistics of random knots obtained by numerical 
simulations with algebraic invariants of knots. In 
particular, we consider how the knotting probability 
should depend on the length TV of ring polymers. 

Let us first consider the knotting probability of 
the trivial knot. We recall that it is the probability 
of forming a trivial knot in a given random polygon 
(or self-avoiding polygon) with N nodes. We denote 
it by Pq {N, r) , and call it the unknotting probabil- 
ity or knotting probability of the trivial knot (or 
unknot). Here, the symbol K — Q denotes that the 
knot K is trivial. We also note that the probability 
Pknotted{N,r) of a polygon being knotted is given 
by the unknotting probability Pq (-^, ''') as in the 
following: Pknotted{N , r) ^ I - Po{N,r). Here we 
recall that if a polygon is equivalent to a nontrivial 
knot, then we regard it as knotted. 

The unknotting probability has been evaluated 
for different models and methods: some models of 
the closedrandom polygons |§, @, js) , the hedge-hog 
method 



13 



For 



and the rod-bead model ||1C 
some models of random polygons or self-avoiding 
polygons, it has been found in Rcf. ||, [l^ that the 
unknotting probability has a decreasing exponential 
dependence on the number N of nodes: 



Po{N)^CoeM~N/No). 



(3.1) 



Here, the two paramete rs N q and Cg are to be de- 
termined so that eq. (3J^) gives the best fitting 
curve to the numerical estimates of the unknotting 
probabilities for different numbers N. We call the 
parameter A^o the characteristic length of random 
knotting of unknot. 

Let us consider the case of nontrivial knots. We 
recall that the knotting probability of a nontriv- 
ial knot K is given by the probability Pk{N) of 
observing the knot if in a given random polygon 
(or self-avoiding polygon) with N nodes. Knotting 
probabilities of several nontrivial knots have been 
evaluated for some different models of random poly- 
gons or self- avoiding polygons. |ll|, |lj, |l5|, |l^ 
Through the simulations with the Vassiliev-type 
knot invariants jl^, |T^, it has been found that, for 
the Gaussian random polygons and the rod-bead 
self-avoiding polygons, the probabilities of nontriv- 
ial knots versus the number N of nodes are well 
approximated by the fitting curves given by the fol- 
lowing formula: 

Pk{N) = Ck {N/NkT'' eM-N/NK). (3.2) 

Here Ck, Nk and ttik are fitting parameters to 
be determined from the estimates of the knotting 



probabilities obtained numerically. For a given knot 
K , we call nriK the exponent of the knot. We also 
call Nk the characteristic length of the knot K. 

There are some universal properties of the fitting 
parameters Nk and uik of the formula (3.2). For 
any one of the models investigated, it is found that 
the values of the parameters Nk are almost the 
same for any knot K: Nk ^ No- @ |l6) Fur- 
thermore, it is observed that the parameter niK 
of a knot K should be universal for the different 
models, while the fitting parameters Ck and Nk 
are model-dependent. Q Thus, we may consider 
that the parameter ttik should play a similar role 
with the critical exponents of critical phenomena. 
In fact, the formula ( |3.2| ) of knotting probability is 
quite consistent with the asymptotic scaling behav- 
iors of the number W{N) of all the configurations 
of self-avoiding polygons: 



W{N) ~ C^i'^N'^ 



for /V > 1 



(3.3) 



Here /i is the growth constant and a corresponds 
to the critical exponent related to the energy of the 
n-vector model through the limit of sending n to 0. 

We can also expect a similar asymptotic expan- 
sion for the number Wk{N) of a knot K. Thus, we 
may call eq. (B^) a scaling formula of the knotting 
probability for the knot type K. 

Let us now consider effects of excluded volume 
on the knotting probability. Through the nu- 
merical simulations of the cylindrical self-avoiding 
polygons, it has been found that the excluded- 
volume parameter such as the cylinder radius plays 
a central role in the probability of knot forma- 
tion. ^, |ll|. Here, we recall that the probability 
Pknotted{N, r) of a Cylindrical self-avoiding polygon 
being a nontrivial knot depends on the number N 
of polygonal nodes, the cylindrical radius r, and the 
Kuhn statistical length h. It was discussed that the 
probability Pknotted{,N , r) should depend on the ra- 
dius r through the exponential factor exp(— 7r/6) 
such as shown in eq. ( |l.2| ), where 7 is the propor- 
tionality constant determined from the simulations. 

am. 



The empirical formula (L2) has been fundamen- 
tal in the study of random knotting of small circular 
DNAs. [|l^ [2^ However, it is not valid when N is 
large enough, as we shall see in §4.3. Here, we recall 
that N is the number of nodes of cylindrical self- 
avoiding polygons. Furthermore, we show in §4.4 
that our generalized scaling-type formula for the 
knotting probability of the cylindrical self-avoiding 
polygons generalizes the formula (|l.2|). 

3.2 Generalized scaling formula for 
knotting probability of ring poly- 
mers 

Let us discuss the knotting probability Pk{N, r) of 
the model of cylindrical self-avoiding polygons con- 
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structed by the dimerization method that depends 
not only on the number TV of nodes but also on the 
cylinder radius r. Generalizing the scaling formula 
( |3.2[ ) which explains only the A^-dependence, we in- 
troduce the following formula for the probability 
PK{N,r): 

PK{N,r) = CK{r) (N/NKir))""'^'^ exp(-iV/iVK(r)) 

(3.4) 

Here Ci<-(r), mK{r), Nxir) are fitting parameters. 
As we shall see later in this section, eq. (^^) gives 
good fitting curves to numerical estimates of the 
knotting probability of the cylindrical self-avoiding 
polygons constructed by the dimerization method. 

The formula (3.4) also generalizes another for- 
mula of the knotting probability. In a previous 
paper |24 , we have discussed how the probability 



of unknot should depend on the cylinder radius r 
for the self-avoiding polygons consisting of cylin- 
drical segments with radius r. We note that the 
cylindrical self-avoiding polygons in Ref. ||2^ are 
constructed by the same method discussed in §2. 
Through numerical simulation, we have found that 
the unknotting probability denoted by Pq (iV, r) can 
be well described by the following: 



Po{N,r) = Co{r) exp {~N/No{r)) , 



(3.5) 



where Co(r), A^o('') are fitting parameters. We note 
that this formula ( |3.5| ) corresponds to the special 
case of the general formula (3.4) when K — Q and 
(r) = 0. We note that the radius-dependence of 
the characteristic length Nq (r) of unknot can be ap- 
proximated by an exponential function of the cylin- 
der radius r p^ : 

A^o(0 «^o(0)exp(ar). (3.6) 

We also note that the exponential dependence is 



favorable to the standard theory of polymers. |29 



3.3 



Numerical estimates of knot- 
ting probability for the model of 
cylindrical self-avoiding polygons 



Let us discuss the numerical data of the knotting 
probability obtained in our numerical simulations. 
For a given number N of nodes, we construct M 
cylindrical self-avoiding polygons with cylinder ra- 
dius r by the dimerization method. Here M is 
given by M = 10*, the polygonal nodes N from 
50 to 1000 by 50 or 100, and the radius r from 0.0 
to 0.007 by 0.001. (For details, see Figure Cap- 
tions.) We have obtained numerical estimates of 
the knotting probabilities for the trivial knot, four 
prime knots (3i,4i,5i,52), and two composite knots 
(3ipi,3ipipi). We note that we shall also dis- 
cuss the knotting probabilities of small sized self- 
avoiding polygons such as = 20 or 30 by the 
cylindrical ring-dimerization method. In §4, the 
connection to the hedge-hog method is studied. 



In Fig. 4, the estimates of the knotting probabil- 
ity Po (iV, r) for the trivial knot are plotted against 
the number N of nodes with r=0.001, 0.003, 0.005. 
The lines in Fig. 4 are theoretical curves given by 
the formula ( |3.5D . We confirm the exponential de- 
cay of Po{N,r) also for the cylindrical self- avoiding 
polygons with large numbers of nodes which are 
constructed by the dimerization method. (See also 
Ref.||.) 

In Fig. 5, the numerical values of the knotting 
probabilities Pk {N, r) are plotted against the step 
number N for the prime knots 3i, 4i, 5i. We see 
from Fig. 5 that the majority of nontrivial prime 
knots are given by the trefoil knots (3i). We also 
observe that the peak values of the knotting prob- 
abilities decrease as the cylinder radius r increases. 
There is also a tendency for the prime knots that 
knots with larger crossing number have less peak 
values. 

In Fig. 6, we show our estimates of the knotting 
probabilities of knots 5i and 52 versus the number 
N of nodes. We see that the knotting probability of 
the knot 52 is always larger than that of the knot 5i. 
However, for both the two knots the knotting prob- 
abilities have their maximal values at almost the 
same number N of nodes. Here we note that the 
error bars of Fig. 6 correspond to half of the values 
of the standard deviations of the knotting proba- 
bility: In Figs. 4, 5 and 7, the error bars present 
the standard deviations of the knotting probability, 
while in Fig. 6 they show half of the values since 
they are too large to be depicted. We discuss how 
we have evaluated the "standard deviations" later 
in §3.3. 

In Fig. 7, the knotting probabilities of two com- 
posite knots 3ipi and 3ipipi, are plotted against 
the number N of nodes for the values of the cylinder 
radius r given by r=0.001, 0.003, 0.005. The knot- 
ting probabilities of the two composite knots have 
their maxima at different values of the number N 
of nodes, which are also different from that of the 
prime knots. 

In Figs. 5-7 of the knotting probabilities, all the 
fitting curves are given by the formula (^^). From 
the Figures, we see that it gives good fitting curves 
to the graphs of the knotting probability Pk {N, r) 
versus the number N of nodes for the cylindrical 
self-avoiding polygons with the different values of 
the radius r. The least-square estimates of the pa- 
rameters mK{r), Cxir), and Nxir) for the fitting 
curves are listed in Table || together with their 
values. The errors for the best estimates of the pa- 
rameters given in Table || correspond to 68.3 % con- 
fidence intervals, which are equivalent to the stan- 
dard deviations. 

Let us explain the method for evaluating the er- 
rors of the estimates of the knotting probabilities: 
with the error values, we have determined the fitting 
curves shown in Figs. 4-7. Throughout the paper. 
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we estimate the error of the knotting probability 
PK{N,r) by the following method: taking the sum 
of the contributions from the statistical fluctuation 
of the number Mk and that of statistical weight 
appearing in the ring-concatenating procedure (i.e., 
the ring-dimerization process in §2.3), multiplying 
the sum by the factor of 2, and then we regard the 
result as the error corresponding to the standard 
deviation of the knotting probability Pk {N, r) . It 
seems that this method of evaluating errors might 
give larger values of errors than other methods. 
However, we have employed it in order to compen- 
sate some possibly neglected contributions to the 
errors arising from the chain-dimerization process. 
It would be not easy to estimate the possible errors 
in the dimerization process, systematically. In fact, 
we do not take into account any possible statisti- 
cal fluctuation in the dimerization process for con- 
structing cylindrical self-avoiding chains. It could 
be as large as that of the statistical weight in the 
ring-concatenating procedure. 

We make a comment on the dimerization pro- 
cess. Giving one initial number to the generator 
of pseudo-random numbers, we construct M — 10* 
cylindrical self-avoiding polygons with number N of 
nodes and radius r. This method makes our simula- 
tion simpler and more practical than giving several 
different initial numbers to the generator. However, 
it might give slightly larger values of errors. This 
consideration has been reflected in our method for 
estimating errors discussed in the last paragraph. 

3.4 Some novel properties of fitting 
parameters 

Let us discuss the best estimates of the fitting pa- 
rameters NK{r), mK{r) and CK{r) listed in Table 
^. From the list, we can derive a number of impor- 
tant observations for the fitting parameters. 

First, we consider the values of the parameter 
CK{r). From Table |, we see that the values of 
Ck {f) are almost independent of the cylinder radius 
r, within the error bars. Thus we may have the 
conjecture that the parameter Ci<-(r) should be a 
constant with respect to the radius r for any knot 
K: 

CK{r) ~ Ck(0). (3.7) 

The parameter (7^(0) of a knot K must reflect 
at least some kind of complexity of the knot. In 
fact, the probability of observing the knot type K 
is small when the value of Cx(0) is small. Further- 
more, as far as the prime knots are concerned, there 
seems to be a tendency that a knot with a smaller 
value of the knotting probability is likely to be more 
complex. Furthermore, the value of the parameter 
Cif (0) is not directly related to the crossing num- 
ber of the knot K . For instance, the value of 6*5^ (0) 
is smaller than that of C52 (0) although both of the 
knots 5i and 82 have the same crossing number. 



Second, we consider the characteristic length 
-/Vif(r). Let us discuss how the parameter Nxir) 
should depend on the cylinder radius r for a given 
knot type K. From Fig. 8, we see that Nxir) 
can be approximated by an exponential function of 
the chain thickness r : Nx{r) = Nfc{0) exp{aKr), 
at least for the trivial knot K = and the trefoil 
knot K = 3i. Considering the poor statistics for the 
case of the nontrivial knots, we may conjecture that 
Nk{t) should depend on r exponentially for any 
knot K. Let us next discuss the knot-dependence 
of the parameter Nk{t) with the radius r fixed. 
From Fig. 9, we see that the values of the parame- 
ters Nxir) for the different knots are almost given 
by the same value with respect to the errors. 

Combining the results of the radius- and knot- 
dependence of the parameter Nxir), we have the 
conjecture that the parameter N^ir) for a knot K 
should be equal to that of the trivial knot and also 
that it should be given by an exponential function 
of the cylinder radius r: Nxir) = Nq{0) exp(ao'')- 

In fact, in the previous paper |Q, we calculated 
the probability of observing the trivial knot from 
the values of the cylinder radius r from 0.0 to 0.1 
by 0.01, and discussed how the parameter iVo(r) 
should depend on the cylinder radius r. Then, 
we found that the characteristic length No{r) is 
roughly approximated by an increasing exponen- 
tial function of r as No{r) = A'^o(O) exp(ar), with 
iVo(O) = 292 ± 5, and a = 43.5 ± 0.6. We note 
that x^=42 while the number of data points is ten 
with the two fitting parameters. In order to improve 
the value, we have also considered another fit- 
ting formula iVo(r) = A^o(O) exp(/3r''), where there 
are three parameters to fit iVo(O), /3, and v. The 
best estimates of the fitting parameters are given 
by A^o(0)=271±6, /3 = 29 ± 2, 1/ = 0.85 ± 0.02, and 
the value is given by 2.2 . 

The results of the previous paper are consis- 
tent with that of the trivial knot in this paper. We 
note that the range of the cylinder radius r in the 
present paper is much narrower than that of Ref. 
p^ . Here, we consider the values of the cylinder 
radius from 0.001 to 0.007 by 0.001. Even the sec- 
ond fitting function: iVo(r) = A^o(O) exp(/3r'^) can 
also be approximated by an exponential function of 
r. However, the connection between the formulas 
(3.4) and (|l.2|)(or (4.1)) can be shown more clearly 



by the first fitting formula than by the second one 
as we shall show in §4. Thus, we employ the first 
one in the paper. 

Third, we discuss the exponent m^f (r) of a knot 
K. From Table |, we see that the values of the pa- 
rameter rriKir) should be independent of the cylin- 
der radius r for any knot K . This property is con- 
sistent with the conjecture that the exponent uik 
of a knot K should be universal, which was also in- 
vestigated for the rod-bead model of self-avoiding 
polygons and the Gaussian model of random poly- 
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gons. ||l6[ Furthermore, within the error bars, the 
values of the exponents can be roughly approxi- 
mated as follows: TO/f(r) ~ for the trivial knot 
{K = 0); uiKir) ~ 1 for the four prime knots 3i, 4i, 
5i and 62; mx('") — 2 for the composite knot 3itJ3i; 

(r) ~ 3 for the composite knot 3ipipi. It is 
interesting to note that the roughly approximated 
values of (r) are consistent with the results not 
only of the off-lattice models such as the Gaussian 
model Q and the rod-bead model Q, but also 
of the lattice model of self-avoiding polygons. 
The roughly approximated values are also consis- 
tent with the additivity of the exponent itik that 
the exponent of a composite knot should be given by 
the sum of the exponents of the constituent prime 
knots: m{Ki^K2) = m{Ki) + m{K2) and so on, 
which was first observed for the Gaussian model of 
random polygons. Q 

Let us discuss an interesting consequence of the 
conjectures on the fitting parameters given in the 
above. We first note that among the three fitting 
parameters, only the characteristic length Nxir) 
should depend on the cylinder radius r; the other 
two parameters should be independent of it. We 
recall also the conjecture that Nxir) — No{r) for 
any knot K. Thus, if we assume that the value of 
the exponent niK can be approximated by the same 
value for any prime knot K, then we have the fol- 
lowing relation between the knotting probabilities 
of any two prime knots Ki and K2: 



Pk, {N, r)/PK, {N, r) ~ Ck, /Ck, ■ 



(3.8) 



Thus, the ratio of the knotting probabilities of any 
two prime knots Ki and K2 should be roughly given 
by that of the parameters Cki and ■ 

4 Topological Entropy and 
Random Knotting of Circu- 
lar DNAs 

4.1 Probability of random knotting 
of circular DNAs 

In this section, we discuss the connection of the 
knotting probability to some experiments of DNA 
random knots |20[ |. 

Let us consider DNA molecules in an electrolyte 
solution. The DNAs are polyelectrolytes with neg- 
ative charges, and the DNA chains are surrounded 
by some clouds of counter ions. Partially due to the 
electrostatic repulsion among the chains, they are 
considered as stiff chains; in fact, the Kuhn statisti- 
cal length of DNA chains is given by rather a large 
value, about 50 nm Q. 

In order to study thermodynamic properties of 
DNA chains we can simulate the DNA chains by 
some configurations of wormlike chains with the 



efi^ective diameter corresponding to the screening 
length [^1], ^ . For any wormlike chain, the persis- 
tent length (i.e., the Kuhn length) is fundamental. 
The length of a wormlike chain should be expressed 
in terms of the Kuhn statistical unit. Thus, we can 
replace wormlike chains by such self-avoiding walks 
that have cylindrical segments with the persistent 
length. Let us now consider wormlike rings which 
are given by rings made of wormlike chains. Then, 
we can also approximate wormlike rings by such 
self-avoiding polygons consisting of cylindrical seg- 
ments with some effective diameter whose length is 
given by the Kuhn length [ pT| . 

Let us discuss the algorithm for generating cylin- 
drical self-avoiding polygons which we call the 
hedge-hog method of ring polymers. The al- 

gorithm is given by the following: we generate a 
set of vectors of unit length with a common ori- 
gin (a "hedge-hog"), and applying the Monte-Carlo 
procedure to the set, we derive a random sequence 
of self-avoiding polygons, keeping only those con- 
figurations that have no overlap between any two 
unadjacent cylindrical segments. (See §A.2) 

Through the numerical simulations, the proba- 
bility of being knotted has been evaluated for the 
hedge- hog method |11 . For the data, the radius- 
dependence of the probability Pknotted{N , r) has 
been well approximated by the empirical formula 
( |l.2| ). Furthermore, a different method of wormlike 
polygons has been introduced, which we call the 
MC method with the bending energy (See §A.3) 
]T9[ , and the knotting probabilities of two prime 
knots 3i and 4i are evaluated for the wormlike rings 
with N Kuhn lengths, where N is given from 16 to 
60. For a given nontrivial knot K , let the sym- 
bol Pk [N, r) denote the knotting probability of the 
wormlike rings, which are constructed by the MC 
method with the bending energy with the chain ra- 
dius r consisting of A'' units of the Kuhn length h. 
Then, the data are fitted by the following formula: 



PK{N,r)^PK[N,Q)e^],{- 



(4.1) 



Here, "fK is a constant which is independent of N 
and to be determined by the simulation for each 
nontrivial knot K ]l9[ |. 

The knotting probability of the wormlike rings of 
the MC method with the bending energy is con- 
sistent with that of the hedge-hog method; the 
knotting probability of wormlike rings of the MC 
method with the bending energy consisting of N 
Kuhn units with the cylinder radius r should be 
equivalent to the hedge-hog method with N cylin- 
drical segments with radius r. In fact, it was dis- 
cussed in Ref. ||l^ that the knotting probability of 
a particular knot for the wormlike ring can be de- 
termined by the number N of Kuhn statistical units 
and the chain radius r. In general, for any model of 
wormlike rings, the Kuhn statistical length b is not 



7 



necessarily equivalent to the length of one segment 
of the rings; it can be much longer than the length 
of the segments, such as the case of the MC method 
with the bending energy. 



The empirical formulae (1.2) and (4.1) give the 



connection of the knotting probability to experi- 
ments of DNA random knots. V.V. Rybencov et.al. 
and S. Y. Shaw and J.C. Wang measured the frac- 
tions of knotted species generated by the method 
of random ring closure of DNAs through the elec- 
trophoretic separation method. ||l^, ^ The effec- 
tive diameter of the DNA double helix was evalu- 
ated, by comparing the knotted fractions observed 
in the experiments with the theoretical estimates 
of the knotting probability given by the computer 
simulations ||ll[ of the hedge-hog method. [^9, ^ 

4.2 Consistency of knotting prob- 
abilities for dimerization and 
hedge-hog methods 

We now discuss that the hedge-hog and dimer- 
ization methods for constructing cylindrical self- 
avoiding polygons are consistent as far as their 
knotting probabilities are concerned. Let us con- 
sider our numerical estimates of the probability 
Pknotted {N, r) of being knotted for the hedge- hog 
and dimerization methods, which are obtained in 
our simulations. Here, we set & = 1 for the hedge- 
hog method without any loss of generality. 

For the hedge-hog method, we have evaluated the 
probability Pknotted{N, r) for the cases of A^=20 and 
30. For the dimerization method, we have evaluated 
it for the cases of iV=21 and 31. For each of the 
methods, the values of the diameter 2r are given by 
the values from 0.0 to 0.05 by 0.01. 

Our numerical estimates of the probability 
Pknotted{N, r) of the hedge- hog method gives almost 
the same with that of Ref. |l^. Furthermore, our 
estimates of the probabihty Pknotted{N , r) of the 
hedge-hog method is also consistent with that of 
the dimerization method. We have thus confirmed 
that the data of the probability PknoUed{N,r) of 
the cylindrical self-avoiding polygons of the dimer- 
ization method satisfy the empirical formula (1.2) 
for the thickness-dependence of the knotting prob- 
ability, for the small numbers N of nodes. 



4.3 



Limited vahdity of the empirical 
formula for the knotting proba- 
bility of DNAs 



We recall that the empirical formula (4.1) describ- 
ing the thickness-dependence of the knotting prob- 
ability has been fundamental in the study of the 
random knots of circular DNAs. However, we shall 
show that it is not valid when the number N of 
nodes is large enough. 



Let us consider the iV-dependence of the ratio: 
PK{N,r)/PK{N,0) of the knotting probability for 
a given nontrivial knot K. If the empirical formula 
(4.1) should be valid, then the ratio would be given 
by the following: 

Pk{N, r)/PK{N, 0) = eM—fKr/b) . (4.2) 

We note that it should be constant with respect to 
N. 

In Fig. 11, the data of Pk{N ,r) / Pk{N,Q) for 
the region of r from r=0.001 to r=0.007 by 0.001 
are plotted against the number N of nodes. Then, 
we see that the ratio is not constant but increases 
with respect to N . Thus, the formula ([4.l[ ) does not 
hold when the number N of nodes is large enough 
such as TV > 100. 

We make a comment on the method for estimat- 
ing the errors shown in Fig. 11. We note that the 
error bars denote the standard deviations in Fig. 
11. The variance ct^ of the ratio: P{N ,r) / P{N ,Q) 
is given by the following formula: 



1 



P{N,QY 



fP{N, r) 
VWO) 



(4.3) 



Here fii denotes the variance of P{N, r) and a2 that 
of P(iV,0). 

4.4 Reduction of the scaling formula 
of knotting probability in the 
case of small 

Let us explicitly discuss how the generalized scal- 
ing formula (^^) is related to the empirical formula 
(4.1). In fact, this gives a kind of "extrapolation" 
of the scaling formula into a region of small N . 
The scaling formula should be valid when N is very 
large, since it is nothing but an asymptotic expan- 
sion with respect to N. On the other hand, the 



empirical formulae (1.2) and (4.1) should be valid 
only for some small number N of nodes. 

There are three points in our discussion. First, 
we recall that the parameter Ck (r) should be inde- 
pendent of the cylinder radius r: Ck{t) — (7^(0) . 
Then, from the formula ( |3.4| ), we have the following 
expression for the ratio of Pk{N, r) to Pk{N, 0): 



PK{N,r) ^ f NKjO) 
Pk{N,0) \NK{r) 



exp 



N 



N 



,7Vo(0) iVo(r) 
(4.4) 

Second, we can neglect the exponential factor in 
the right-hand-side of (4.4), when N is small such 
as < 30. The characteristic length N^ir) can 
be roughly evaluated as Nx{r) > 300 when r is 
given by some value from 0.001 to 0.01. The value is 
about ten times larger than that of the number N of 
nodes, when N < 30. Therefore, the exponential ar- 



gument 



N 



N 



should be rather small. Third, 
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we recall that the characteristic length No{r) can 
be given by A^o('') — -^o(O) exp(ar), as discussed in 
§3.4 (see also in the previous paper: Q). Further- 
more, we may assume the conjecture that the char- 
acteristic lengths do not depend on any knot type 
K: Nj^{r) ~ A'o(r). Combining the three points 
given in the above, we see that the main contribu- 



tion to the ratio of eq. (i.4) is given by 



PKjN, r) 
Pk{N,0) 



• exp {—arriKr) 



(4.5) 



We note that the factor (iVK(0)/iVK(r))'"^- should 
correspond to the exponential factor exp(— 7A'r/6) 
in the empirical formulae (1.2) and (4.1). Explicitly, 
we have = OirriK when h — 1. Thus it is shown 
that the empirical formula (1^) can be derived from 
the generalized scaling formula (3^). We recall that 
the constant is assumed to be independent of N 
in Ref. 

We remark that the knotting probability 
Pk {N, r) can be described by the generalized scal- 
ing formula (3^) from small N to large N , as shown 
in §3.3, explicitly. Furthermore, if we assume the 
properties Ci^(r) ~ Ck(0), Nxir) ~ Noir) and 
No{r) ~ A^o(O) exp(ar), then we have the following 
expression 



that for any knot investigated, the knotting proba- 
bility of the model of cylindrical self-avoiding poly- 
gons is described by the generalized scaling formula 
(3.4) as a function of the number N of nodes and 
the cylinder radius r. 

From the best estimates of the fitting parameters, 
we have observed several important properties of 
the parameters. Based on these properties, we can 
derive a conjecture of the best formula of the knot- 
ting probability. For self-avoiding polygons with 
N cylindrical segments with the radius r of unit 
length, the best formula of the knotting probability 
for knot K is given by the following: 

PK{N,r) = Ck [N/N.ir))'^"' exp(-Ar/7V,(r)) , 

(5.1) 

where the characteristic length Nc{r) is indepen- 
dent of knots and is given by 



iV,(r) =7V,(0)exp(ar) 



(5.2) 



log 



PK{N,r) 
PKiN,0) 



N 



-aniKr- 



NoiO) 



(l-exp(-ar)), 
(4.6) 



which is consistent with the A^-dependence shown 
in Fig. 11, within error bars. 

For the case of the trivial knot, we may also con- 
sider the r-dependence of the second factor of the 



right-hand-side of eq. (4.4), which can be neglected 
for the case of nontrivial knots. Taking into ac- 
count the fact that too ~ 0, we have the following 
approximation: 



We recall that Ck and tjik should be constant with 
respect to the cylinder radius r. We shall investi- 
gate the conjecture more precisely in later publica- 
tions. 

Finally, we make a comment on the range of r. 
In the paper, we have discussed for some nontrivial 
knots the knotting probabilities of the cylindrical 
self-avoiding polygons with radius r given from to 
0.007. On the other hand, in Ref. Q, we have dis- 
cussed the probability of unknot for the wider range 
of r: from to 0.1. From the result of Ref. ||2^ , it is 
suggested that the conjecture (5.1) should be vahd 
also for the larger values of radius r. Thus, it is an 
interesting future problem to check the conjecture 
explicitly through numerical simulations. 
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Appendix A: Algorithms for the Model of 
Ring Polymers 

In Appendix A, we describe the three different 
algorithms for constructing self-avoiding polygons 
consisting of freely jointed hard cylinders. |2j, ^ 
^ All the three algorithms should produce equiva- 
lent sets of ring polymers corresponding to the dis- 
crete worm-like chains, if the numbers of polygons 
constructed are very large. As far as constructing 
large polygons, however, it seems that the dimer- 
ization algorithm should be the most efficient. 

A.l Cylindrical 

ring-dimerization method (the 
dimerization method) 

Let us explicitly discuss the method for generating 
ring polymers of freely jointed hard cylinders, by 
which almost all the polygons in the paper are con- 
structed. The method is based on the algorithm 
of Ref. [^ for the rod-bead model. It consists of 
three parts: (1) generation of basic chains of hard 
cylinders; (2) propagation of linear chains by the 
dimerization algorithm; (3) formation of ring poly- 
mers by connecting two linear chains. We call the 
combined procedure of the parts (1), (2) and (3), 
the cylindrical ring-dimerization method. 

Basic chains consisting of hard cylinders are gen- 
erated by the straightforward Monte Carlo method. 
For an illustration, let us suppose that we make 
basic chains of eight cylindrical segments. Then, 
for any pair of unadjacent cylinders, we check the 
overlapping condition between the cylinders which 
is described in Appendix B. If there is a pair satisfy- 
ing the overlapping condition in a given basic chain, 
then we throw away the whole chain and start from 
the beginning. If we find that there is no overlap- 
ping pair of unadjacent cylinders, then we store the 
chain in a computer disk. Here we recall that we do 
not check the overlapping condition for any pair of 
adjacent cylinders. This corresponds to the freely- 
jointedness of the model. 
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Let us consider how to make chains of sixteen 
cyhnders. We make two pools of M basic chains in 
the disk. Then, we choose randomly basic chains 1 
and 2 from the first and second pools, respectively. 
We join them together by placing the zeroth node 
of the second chain on the eighth node of the first 
chain. If there is no overlapping pair of unadjacent 
cylinders in the new chain of sixteen cylinders, then 
we store it in the computer disk. If there is any 
overlapping unadjacent pair, then we throw it away. 
We thus construct M chains of sixteen cylinders in 
the disk. The same procedure is used for creating 
chains of , say, 32 cylinders, 64 cylinders, etc. We 
note that chains with any number of cylinders can 
be constructed by the dimerization scheme shown 
in Fig. 3. 

Let us describe the procedure of making ring 
polymers of 2N + 1 cylinders from linear chains of 
TV cyhnders. We consider two pools Si and 52 of 
M chains of freely jointed N hard cylinders. We 
pick up randomly chain 1 in the pool Si . Then, we 
can choose chain 2 randomly in the pool S'2. In or- 
der to make the process more efficient, however, we 
choose chain 2 as follows. Let the symbols hi and 
ft,2 denote the end-to-end distances of chains 1 and 
2, respectively. Here we recall that the length of 
cylindrical segments is given by 1. Then, we choose 
chain 2 randomly from the group of chains in S'2 
which satisfy the condition: 



Here, the angle 9 is determined by the values hi 
and /i2 as follows 



\hi - 11 < h2 < hi + 1. 



(A.l) 



The sampling bias induced by this operation is cor- 
rected with the statistical weight m/M, where m is 
the number of chains in S'2 which satisfy the condi- 
tion (A.l). We recall that M is the total number of 
linear chains in S2. 

Let us put a new cylindrical segment between the 
two ends of the chains 1 and 2. Then, we check 
whether there exists any segment in the two chains 
which overlaps with the new segment, except for the 
two neighboring ones. We also check the overlap- 
ping condition for every unadjacent pair of cylinder 
in which are cylinder is in chain 1 and another in 
chain 2. If there is no overlap, then we consider that 
the selected pair of chains makes a perfect ring with 
2A^ -I- 1 cylindrical segments. We store the ring in 
the hard disk. However, the probability of forming 
a ring depends on the values of /ii and /i2 . Let us de- 
note by 6 the angle between the end-to-end vectors 
of the two chains 1 and 2. Then, the probability is 
proportional to 27rsin6'. We note that 27r/i2sin0 is 
the arclength of the circle on which the end-point of 
the chain 2 can be placed. The probability of form- 
ing a ring should be proportional to the arclength of 
the circle divided by ft,2. Thus, the total statistical 
weight of the ring-dimcrization procedure is given 

by 

W^—sinO (A.2) 



cos 9 



hi 



hi 



1 



2/li/l2 



(A.3) 



We remark that all the expectation values of some 
quantities of the model such as the knotting prob- 
ability etc., should be calculated by taking the 
weighted averages with respect to the statistical 
weight of eq. ( |A.2| ). 

A.2 Monte-Carlo method with ran- 
dom hedgehog (the hedge-hog 
method) 

Let us describe the algorithm of Ref. [jll] , which we 
call the Monte-Carlo method with random hedge- 
hog, or the hedgehog method, for short. First, a 
set of n (even number) vectors of unit length 
with a common origin (a "hedgehog") is generated 
as follows: In the hedgehog, for those vectors that 
have odd i values, their directions are chosen ran- 
domly and independently; for even i values, we set 
Ci — — e,j_i. Thus, we have — ^- Then, to 

exclude pair correlations between the vectors, for 
each Ci, we choose a different vector randomly 
from the hedgehog and the pair of vectors and 
Cj are rotated at a random angle around the bi- 
sectrix of the angle between them. This operation 
does not change the sum -I- Cj , and the sum of 
all the vectors in the hedgehog remains zero, conse- 
quently. The process is repeated for each vector of 
the random hedgehog, many times. As a result, a 
"random hedgehog" is obtained. Then, the vectors 
of the random hedgehog are put in an arbitrary or- 
der, and we have the chain. The resulting chain is 
automatically closed. 

In order to take into account the excluded volume 
effect, chain segments are modeled as hard cylinders 
with radius r. We generate a large number of closed 
chains by the random hedgehog method, and then 
retained only those of them for which the minimum 
distance between any two unadjacent segments is 
larger than the segment diameter 2r. 

A.3 Monte-Carlo method with the 
bending energy 

We consider the algorithm of Ref. ||l^, which we 
call the Monte-Carlo method with the bending en- 
ergy. A closed chain of kn segments of rigid im- 
penetrable cylinders of equal length and radius r 
is constructed. Here k elementary segments corre- 
spond to the Kuhn statistical length, and the closed 
chain consists of n Kuhn segments. The conforma- 
tional sets are obtained by successive deformations 
of a starting conformation in accordance with the 
Metropolis-Monte Carlo procedure. The deforma- 
tion is rotation of a subchain containing an arbi- 
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trary number of adjacent segments by a randomly 
chosen angle, 0, around the straight line connecting 
the vertices bounding the subchain. The value of (p 
is uniformly distributed over an interval {—(po to 
(j)o), where the interval is chosen so that about half 
of the moves are accepted. It depends on the energy 
whether a trial conformation generated is accepted 
or not. The energy of the chain, E, is calculated as 

kn 

E = RTa^e1 (A.4) 

i=l 

where summation is done over all the joints between 
the elementary segments, R is the gas constant, T 
is the absolute temperature, 9i is the angular dis- 
placement of the zth segment relative to segment 
z — 1, and a is the bending rigidity constant. The 
bending rigidity constant is chosen so that exactly 
k elementary segments correspond to the Kuhn sta- 
tistical segment length. For each set of values of n 
and r, a large number of conformations are gener- 
ated. 

Appendix B: Overlapping Condition of the 
Cylinder Model 

We discuss the condition when a given pair of 
cylinder segments with radius r has an "overlap" , 
explicitly. We first recall that the central line seg- 
ment of a cylinder is defined by the line segment 
between the centers of the upper and lower disks 
of the cylinder. We now define the condition of an 
"overlap" as follows: given two cylinder segments 
are said to have an "overlap" if and only if the 
distance between their central line segments is less 
than the cylinder diameter 2r. 

Let us formulate the algorithm for the overlap- 
ping condition. We consider a pair of two cylinders 
with radius r and unit length in three dimensions. 
We may assume that the end points of their central 
line segments are given by b and a + b, d and c + d, 
respectively. Here, the vectors a and care unit vec- 
tors. Then, any point on the central line segments 
can be expressed by Xg or Xt given in the following: 



Xs = sa + b Xt = tc + 



(B.l) 



where s and t are real parameters satisfying < 
s,t < 1. We also define the angle parameter 9 by 
the relation: a-c = cos 6, where we take the branch: 

< 61 < TT. 

Let us consider the two infinite lines extending 

the two central line scgmcnits for Xg and Xt, re- 
spectively. Then, the distance between the two in- 
finite lines is given by the minimum of the quantity: 
D —\ Xs — Xt I . Denoting the minimum distance 
by Drain, we have 

DLn = - tin + 2s,„t™ COS + 7 , (B.2) 

where the parameters and tm are given by 



Here a, /?, and 7 have been defined by a = —a • 
(6 — d), P=c - {b — d), and 7=(6 — d)"^, respectively. 
We note that Sm and tm are such values of the pa- 
rameters s and t that give the minimum Dmin for 
D=\X,-Xt |. 

Let us now explain the algorithm. First we cal- 
culate the value of cos 9. Here, we suppose that 
cos^ < 1, and we shall consider the case when 
cos^ 9=1, later. If Dmin is larger than the diam- 
eter 2r, then there is no overlap between the two 
cylinder segments. If Dmin is smaller than 2r, we 
check the values of Sm and tm- If 0< Sm <1 and 
0< tm <1, there is an overlap; otherwise, we choose 
the value dl^in, which will be given shortly, and we 
calculate £)^„„ = d^j„ -|- £)^i„. If Dsum is larger 
than 2r, then there is no overlap; otherwise, there 
is an overlap. 

Let us now formulate the value rf^j„, explicitly. 
We first define for j = 0, 1, 2, 3, in the following: 



ZO = Sm- tm COS ( 

Zl = tm - Sm COS ( 

Z2 = Zq + COS 9 

Z3 — Zl + COS 9 



(B.4) 



Pi = < sisin^ 



Then, we define for j = 0, 1, 2, 3 as follows: 

+ t^m sin^ 9 for zo < 

Pq = <^ t^ siv? 9 for < zo < 1 

(1 - Zof + 4 sin^ 9 for zq > 1 

zl + sin^ 9 for zi < 

for < Zl < 1 
(1 - zi)2 + s^^ sin^ 6* for zi > 1 

zl + {l- tmf sin^ 9 for Z2 < 

P2 = \ {I- tmf sin^ 9 for < Z2 < 

(1 - Z2f + (1 - tmf sin^ 9 for Z2 > 1 

z| + (1 - Smf sin^ 9 for Z3 < 

P3 = <; (1 - Smf sin2 9 for < Z3 @ 

{1 - zaf + {1 - Smf sin'' 9 for Z3 > 1 

Then, we define dmm by the minimum of pj for 
J = 0,1, 2, 3. 

Finally, we discuss the case when 9 = or tt. 
When 6* = 0, we define -D^j„ by 



1 -h 2a -h 7 for a < -1 



7 - 



for -1 < a < 2 (B.6) 



1 



sin 



-{a+P COS 9) , 



1 



sin 



-{acos9+P) 
(B.3) 



1 - 2a + 7 for a > 2 
When 9=n, we define it by 

{7 for a < 

7-a2 for 0<a<2 (B.7) 

4 - 4a + 7 for a > 2 

If Dmin is greater than 2r, we have no overlap; oth- 
erwise there is an overlap. 
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Figure Captions 



Fig. 1: Trivial knot, prime knot 3i (trefoil knot), 
and composite knot 3i#3i. Knot 3i#3i is given by 
a product of two prime knots 3i. 



Fig. 2: Polygonal knot equivalent to trefoil knot 3i. 

Fig. 3: Dimerization scheme for constructing self- 
avoiding walks of length 100; chains of lengths 12 
and 13 are generated by a direct methods; then, 
chains of 25, 50 and 100 are given by the dimeriza- 
tion method, systematically. 

Fig. 4: Unknotting probability PQ{N,r) versus 

number A'^ of polygonal nodes of the cylindrical self- 
avoiding polygons constructed by the dimerization 
method. Numerical estimates of Po{N, r) are shown 
for the following three values of r by black circles, 
black triangles and black diamonds, respectively: 
(a) r=0.001, 0.003 and 0.005; (b) r=0.002, 0.004, 
and 0.006. Error bars denote their standard devia- 
tions. Number N of nodes are given by 51, 151 and 
lOOj + 1 fori = 1,2,..., 10. 

Fig. 5: Knotting probability PK{N,r) of the cylin- 
drical self-avoiding polygons constructed by the 
dimerization method versus number iV of polygo- 
nal nodes for nontrivial prime knots. Numerical 
estimates of PK{N,r) for K=3i, 4i and 5i are 
shown by black circles, black triangles and black 
diamonds, respectively, for the following values of 
r: (a) r=0.001; (b) r=0.003; (c) r=0.005. Error 
bars denote their standard deviations. Number N 
of nodes are given by 51, 151 and 100 j + 1 for 
j = 1, • • • , 10. 

Fig. 6: Knotting probability PK{N,r) of the cylin- 
drical self-avoiding polygons constructed by the 
dimerization method versus number TV of polygo- 
nal nodes for two knots with five crossings. Nu- 
merical estimates of PK{N,r] for if=5i and 52 are 
shown by black circles and black triangles, respec- 
tively, for the following values of r: (a) r=0.001; (b) 
r=0.003; (c) r=0.005. Error bars depicted in Fig. 6 
are given simply by the sum of the statistical fluctu- 
ation of the number Mk and that of the statistical 
weight of the ring-dimerization procedure. Number 
N of nodes are given by 51, 151, and 100 j + 1 for 
j = l,2,..-,10. 



drical self-avoiding polygons constructed by the 
dimerization method versus number N of polygo- 
nal nodes for two composite knots. Numerical esti- 
mates of PK{N,r) for K = 3ipi and 3ipitt3i are 
shown by black circles and black triangles, respec- 
tively, for the following values of r: (a) r=0.001; 
(b) r=0.003; (c) r=0.005. Error bars denote their 
standard deviations. Number N of nodes are given 
by 51, 151 and lOOj -|- 1 for j = 1, 2, • • • , 10 . 

Fig. 8: Semi-logarithmic plot of characteris- 
tic length Nxir) versus cylinder radius r of the 
cylindrical self-avoiding polygons constructed by 
the dimerization method. Numerical estimates of 
Nxir) for K=0 and 3i listed in Table 2 are depicted 
by black circles and black triangles, respectively, to- 
gether with their error bars, with the values of r 
from 0.0 to 0.007 by 0.001. All the black triangles 
(knot 3i) are slightly shifted rightward by 0.0001 
for graphical convenience. The straight fitting line 
is determined by the least square method. 

Fig. 9: Semi-logarithmic plot of characteris- 
tic length Nxir) versus cylinder radius r of the 
cylindrical self-avoiding polygons constructed by 
the dimerization method. Numerical estimates of 
Nxir) for K=0, 3i, 4i, and 5i listed in Table 2 are 
depicted by black circles and black triangles, blac;k 
diamonds and black crosses, respectively, for the 
values of r from 0.001 to 0.007 by 0.001, together 
with their error bars. Black triangles (3i), black di- 
amonds (4i) and black crosses (5i) are shifted right- 
ward by 0.0001, 0.0002 and 0.0003, respectively. 

Fig.lO: Probability Pknotted{N,r) of being knot- 
ted versus cylindrical diameter 2r for the hedge- 
hog and dimerization methods. Numerical esti- 
mates of Pknotted{N,r) of the hedge-hog method 
for = 20 and 30 are shown by black circles and 
black triangles, respectively. Numerical estimates 
of Pknotted{N,r) of the dimerization method for A'' 
= 21 and 31 are shown by black diamonds and black 
crosses, respectively. Error bars denote their stan- 
dard deviations. 

Fig. 11: The ratio PK{N,r)/PK{N,0) ver- 
sus number N of polygonal nodes of the cylin- 
drical self-avoiding polygons constructed by the 
dimerization method. Numerical estimates of 
PK{N,r)/PK{N,0) for r = 0.001, 0.005 and 0.007 
are shown by black circles, black triangles and black 
diamonds, respectively, for the following knots: (a) 
trivial knot; (b) trefoil knot 3i. 



Fig. 7: Knotting probability Pk {N, r) of the cylin- 
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Table 1: Values of the determinant of knot | 
Ak{—1) I and the second Vassiliev invariant ^^(K) 



for some simple knots. 


Knot K 


1 Ak(-i)| 


V2{K) 





1 





3i 


3 


-12 


4i 


5 


12 


5i 


5 


-36 


52 


7 


-24 


3itt3i 


9 


24 


3itt3ipi 


27 


-36 



Table 2: Fitting parameters mK{r), Cxir), Nxir) 
to the cylindrical self-avoiding polygons constructed 
by the dimerization method 



trivial 


r=0.001 


r=0.002 


r=0.003 


r=0.004 


r=0.005 


1-0.006 


1=0.007 


Cif(r) 


1.055±0.065 


1.048±0.006 


1.066±0.061 


1.034±0.060 


1.001±0.058 


; 1.034±0.057 1.024±0.055 


NK{r) 


278itl4 


298±15 


315±16 


341±17 


377±20 


388it20 


417±21 


rriKir) 


-0.004±0.034 


-0.002±0.03 


0.0010±0.031 


-0.004±0.030 


-0.019±0.028 0.001±0.027 -0.006±0.025 




3.0 


2.7 


1.2 


2.1 


2.0 


2.2 


1.7 


3i 


r=0.001 


r=0.002 


r=0.003 


r=0.004 


r=0.005 


r=0.006 


r=0.007 




0.625±0.021 


0.638±0.020 


0.655±0.061 


0.669±0.019 


0.675±0.020 


0.687±0.021 


0.678±0.020 


NK{r) 


264±21 


280±22 


297±24 


314±26 


341±30 


375±36 


390±38 


mK{r) 


1.06±0.10 


1.091±0.100 


1.098±0.097 


1.128±0.097 


1.082±0.095 


1.048±0.094 


1.072±0.095 




0.9 


1.7 


0.5 


1.0 


1.9 


0.7 


2.6 


4i 


r=0.001 


r=0.002 


1-0.003 


r=0.004 


r=0.005 


r=0.006 


1=0.007 


Cxir) 


0.124±0.009 


0.125±0.009 


0.125±0.009 


0.128±0.012 


0.125±0.009 


0.126±0.009 


0.131±0.008 


Nuir) 


236±40 


259±47 


286±.56 


348±75 


348±79 


347±79 


334±70 


rriKlr) 


1.269±0.236 


1.210±0.211 


1.176±0.21 


0.978±0.21S 


1.118±0.23r) 


1.118±0.213 


1.247±0.2.36 


X' 


2.0 


2.3 


1.0 


2.0 


i.7 


2.3 


1. 1 


5i 


r=0.001 


r=0.002 


r=0.003 


r=0.004 


r=0.005 


r=0.006 


r=0.007 


CK(r) 


0.035±0.006 


0.041±0.005 


0.040±0.005 


0.037±0.005 


0.034±0.007 


0.036±0.004 


0.035±0.005 


NK{r) 


198±56 


246±72 


311±116 


265±87 


222±62 


326±129 


425±222 


mK(r) 


1.612±0.486 


1.335±0.395 


1.103±0.428 


1.380±0.436 


1.797±0.452 


1.306±0.448 


1.118±0.473 


X^ 


1.6 


3.6 


0.7 


2.1 


2.6 


0.9 


0.8 


52 


r=0.001 


r=0.002 


1=0.003 


r=0.004 


r=0.005 


r=0.006 


1=0.007 




0.068±0.007 


0.072±0.007 


0.068±0.006 


0.064±0.006 


0.060±0.006 


0.062±0.006 


0.061±0.006 


NK{r) 


226it49 


230±51 


283it73 


278±73 


322it96 


330itl01 


290±81 


rriK (r) 


1.401±0.331 


1.404±0.340 


1.234±0.324 


1.377±0.347 


1.317±0.340 


1.296±0.348 


1.418±0.367 


X^ 


1.3 


2.0 


1.2 


1.6 


2.0 


1.3 


1.0 



3iii3i 


r=0.001 


1=0.002 


1=0.003 


1=0.004 


1=0.005 


1=0.006 


1=0.007 


CK(.r) 


0.166±0.026 


0.171±0.029 


0.161±0.032 


0.187±0.032 


0.184±0.034 


0.173±0.038 


0.161±0.040 


NK{r) 


259±36 


273±41 


278±43 


326±57 


318±52 


328±58 


325±59 


rriK (r) 


2.150±0.238 


2.178±0.248 


2.311±0.260 


2.102±0.251 


2.225±0.245 


2.327±0.266 


2.448±0.282 


X^ 


1.9 


0.9 


1.4 


1.9 


2.0 


1.1 


1.6 


3itt3itt3i 


r=0.001 


r=0.002 


r=0.003 


r=0.004 


r=0.005 


r=0.006 


r=0.007 


CKir) 


0.025±0.020 


0.023±0.021 


0.049±0.038 


0.021±0.022 


0.044±0.039 


0.030±0.036 


0.026±0.029 


NK{r) 


249±77 


244±78 


356±152 


260±90 


362±163 


330±161 


318±138 


rriKir) 


3.316±0.689 


3.443±0.741 


2.838±0.672 


3.602±0.783 


2.976±0.719 


3.358±0.888 


3.478±0.804 


x' 


0.2 


1.0 


0.5 


0.7 


0.8 


1.3 


0.9 
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compos It knot: 



(a) Probability of nontrivial knots (r=0.001) 

0.25 I ■ I 




200 400 600 800 1000 1200 



N : Number of nodes 




N : Number of nodes 



(a) Probability of composite knots (r=0.001) 





r : radius of cylinder 




100 ' ' 

0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 



r : radius of cylinder 



Probability of being knotted 



X 



f 



X 



i 









• 


n=20 


0.08 






T 


n=30 




■ i 


♦ 


n=21 


0.06 






X 


n=31 



0.02 0.04 0.06 0.08 



0.1 



0.12 



2r : cylinder diameter 



(a)trivial knot 
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(b)trefiol knot 

T— 1 — I— I — 1—1 — I— I — 1—1 — I— I — 1—1 — I— I — 1—1 — I— I — 1—1 — r 
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